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Abstract 

We report on the design and performance of a near infra-red Mueller 
matrix imaging ellipsometer, and apply the instrument to strain imag- 
ing in near infra-red transparent solids. Particularly, we show that 
the instrument can be used to investigate complex strain domains in 
multi-crystalline silicon wafers. 

1 Introduction 

In this paper the Near Infra-Red (NIR) Ferroelectric Liquid Crystal (FLC) 
based Mueller matrix ellipsometer (MME) design reported recently [I], is 
modified to an imaging set-up, and used to demonstrate the application to 
strain imaging. Neither a NIR imaging MME system nor an FLC based 
imaging MME has to our knowledge been previously reported. 

Several non-scanning imaging MME systems based on Variable Liquid 
Crystal Retarders (VLCR) have recently been reported 013]. These systems 
have so far only been operated in the visible range, and will generally be 
many orders slower than an FLC based system [I]. Using such systems, 
it has been reported that imaging MME is an interesting tool in both bio- 
applications [21 HI G3 [6] , an d strain imaging [7]. Recent non-imaging Mueller 
matrix ellipsometric studies of bio-tissue demonstrate the usefulness of the 
Mueller matrix in combination with polar decomposition techniques [HI 

maim. 

The advantage of the complete Mueller matrix measurement relies on 
the possibility to use polar decomposition techniques [HJ [T2l [T3l HU [15] and 
that the Eigenvalue Calibration Method (ECM) may be used for the calibra- 
tion of the system [T6J [TT] . Furthermore, the Mueller matrix will not suffer 
from unexpected properties initially believed not to be part of the sam- 
ple properties, such as polarization dependent scattering or depolarization, 
different types of diattenuation, and both circular and linear birefringence 
components. 



Several FLC based designs have been proposed, although the first pro- 
posed system similar to the one reported here, appears to be by Gandorfer 
et al. [E]. FLC based MMEs are appealing since they involve no moving 
parts, and supply a highly stable beam. They are both suitable for direct 
imaging applications (hyper-spectral and monochrome), or in conjunction 
with e.g. stripe CCD spectrographs commonly used in spectroscopy. Fur- 
thermore, FLCs are fast, and thus allows for a fast determination of the 
Mueller matrix. The disadvantage of liquid crystals, and in particular the 
current FLCs is the well known degradation upon UV-radiation. Applica- 
tions of liquid crystal technology in MME are so far limited to the visible 
and the infra-red. 

For the growing silicon solar cell industry, it is a major concern to reduce 
the material cost. One approach is to reduce the wafer thickness. A factor 
making the latter difficult, is the internal residual strain, which is often 
induced in the process of casting. To control and verify a successful wafer 
production with lower strain, effective instruments is needed to measure the 
residual strain. In addition to a report of the system instrumentation, we 
demonstrate the application of the FLC based NIR MME imaging system 
to make a map of two dimensional projections of strain fields in multi- 
crystalline silicon wafers. 

2 Experimental 

2.1 Overview of the system and concept 
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Figure 1: Schematic drawing of the Mueller matrix imaging system. The 
system consists of two sources, polarization state generator and analyzer, 
imaging optics and a camera. 
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The NIR-MME imaging setup was designed in order to operate both in 
transmission and reflection mode. The applications demonstrated in this 
paper are only from transmission mode measurements. The system was de- 
signed to operate in the range 700 — 1600nm [2 [19]. A schematic drawing of 
the setup is shown in Figure[IJ In particular, the Polarization State Analyzer 
(PSA) and the Polarization State Generator (PSG) are both composed of a 
NIR dichroic polarizer (P1/P2), a zero order quarter wave-plate at 465nm 
(R1/R4), a half wave-plate at 1020nm (R2/R3), an FLC half wave retarder 
at 510nm (F1/F4) and an FLC half wave retarder at 1020nm (F2/F3). The 
FLCs have fixed retardances and can switch between two orientations of the 
fast axis. A suitable combination of wave-plates and FLCs ensures that the 
PSG generates four sub-optimal Stokes vectors for all wavelengths in the 
design range. These four Stokes vectors define the column vectors of the 
PSG matrix (W). The four analyzer states of the PSA define the rows of 
the analyzer matrix (A). A Mueller matrix measurement is carried out by 
measuring the 16 intensities obtained when switching through all possible 
FLC combinations. The measurement and the calibration process is similar 
to the one proposed for the LCVR systems [20], and the rotating compen- 
sator systems |21j . The 16 measured intensities form a matrix B, which 
correspond to the matrix product of A, the sample Mueller matrix (M) and 
W 

B = WMA. (1) 

M can then be found by multiplying B by the inverse of W and A from 
each side 

M = W _1 BA- 1 . (2) 

It is evident from linear algebra that for equation [2] to be solvable, A and 
W need to be non-singular. Specifically, it has been shown that the error 
in the intensity measurements (B) and the calibration errors of A and W 
are amplified into the errors of a measured Mueller matrix as 



I AMI | / IIAAII MAWM IIABII 

~ rnm + M „ rM + ^w^a UTyll , (3) 



||M|| ~ ||A|| vv ||W|, 

where an d ka are the matrix condition numbers of A and W, defined 
in our work by the L2 norm. 

In order to find W and A, an implementation of the ECM is used. The 
method allows W and A to be calculated without exact knowledge of the 
Mueller matrix of the calibration samples. More details about the FLC 
MME system and its calibration are reported elsewhere fT| [T9j. 

2.1.1 Light sources and detectors 

The system is operated with three different light sources. For spectroscopic, 
and for measurement of weakly scattering samples, a tungsten halogen white 
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light source (150W) is used in combination with a grating monochromator 
or band pass filters. For samples with more absorption or scattering, higher 
intensity was needed and a 980nm diode laser (max 300mW) source or a 
LED array with center wavelength at 1300nm (HOmW) was used. In all 
cases, the light is collimated before entering the PSG. In the case of the 
diode laser, a rotating diffuser was used to reduce speckle. 

The digital camera was an Xenics Xeva camera operating at 15Hz, con- 
sisting of a 14-bit InGaAs FPA detector with 640 x 512 pixels and sensitive 
to the spectral band 0.9 — 1.7/xm. In addition, due to less thermal noise 
and dead pixels, a silicon camera from Hamamatsu was used in combination 
with the 980nm laser. The field of view of the system is 1cm x 1cm with a 
resolution of ~ 13/xm. 

2.2 Analysis of the Mueller matrix 

There are several advantages of complete Mueller matrix ellipsometry com- 
pared to generalized ellipsometry, such as the introduction of the ECM, 
which makes exact modeling of the complex system superfluous. Secondly, 
it is a complete measurement, that includes full polarimetric information 
about the sample, retardance, diattenuation and depolarization. The po- 
larimetric information can in many cases be extracted using polar decom- 
position techniques. Several decomposition techniques have been proposed. 
For imaging purposes the most relevant are the forward [12] , reverse [2] In 
particular, the forward polar decomposition technique, described in details 
elsewhere QUI [HJ [12j [14] , was applied to all measurements presented here. 
The basic principle of the forward decomposition is to assume that the po- 
larizing properties of the measured Mueller matrix (M) is taking place in the 
following order, diattenuation (M^), retardation (M#) and depolarization 
(Ma), which give 

M = M A M R M D . (4) 

It is further convenient to calculate a numeric quantity for the retardance 
(5), the orientation of slow axis the degree of polarization (P), the 

diattenuation (D) and the orientation of diattenuation (Op) [10J. For the 
samples shown in this paper, only the retardance and the orientation of the 
slow axis are explicitly used. 

Due to detector noise, the measured Mueller matrix will always have 
an error. The error may result in a slightly unphysical Mueller matrix. 
By putting appropriate physical constraints on the measured matrices, a 
measure of the matrix physicality can be found. It is further possible to 
calculate the closest physical matrix of an unphysical matrix. In particular, 
the constraints proposed by Cloude [22] was applied in this work. 
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3 System performance, results and discussion 



3.1 Validation of calibration of the NIR MME Imaging Sys- 
tem 

Figure [2] shows the spectral condition number of the system, reported in 
more detail elsewhere pQ, and Figure [3] shows the W image after calibration 
obtained with the 980nm source. The average inverse condition number of 
this matrix was 0.46 for both A and W. 

An estimate of the accuracy of the Mueller matrix imaging system can 
be made by considering a measurement of the well defined Mueller matrix 
of air (the identity matrix) . The mean matrix of a measured Mueller matrix 
image was 



1.000 ±0.00 
0.000 ±0.01 
0.001 ±0.01 
0.001 ±0.01 



-0.004 ±0.01 
1.002 ±0.02 

-0.002 ±0.02 
0.002 ±0.02 



-0.003 ±0.01 
-0.001 ±0.02 
1.006 ±0.03 
0.007 ±0.02 



0.006 ±0.01 
-0.002 ±0.02 
-0.003 ±0.02 

0.998 ± 0.03 



with the standard deviation from the mean value given as the ±. The error 
in the mean matrix is in the third decimal, while the standard deviation from 
the mean is in the second. The variations found over the pixels may have 
many sources, but is most likely due to detector noise. The non-uniformity 
and error is found to be satisfactory, although there is room for further 
improvements. 
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Figure 2: The simulated and measured spectroscopic condition number of 
the system, for both A and W PQ. 
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Figure 3: The W matrix image at 980 nm. The four columns of the matrix 
are the four probing Stokes vectors from the PSG. 



3.2 Polarimetric images and Strain mapping 

Polarimetric imaging can be used to investigate strain fields in transpar- 
ent crystals. An isotropic material which experiences a non-isotropic stress, 
becomes, as described by the photoelastic theory [2311211, anisotropic (bire- 
fringent) due to the induced internal strain. Polarimetry is very sensitive to 
changes in the retardance 

6 = An • d = (n s i ow - n fast ) • d, (5) 

which is proportional to the birefringence (difference in refractive index from 
slow to fast axis, An) and the thickness of the sample (d). The refractive in- 
dex will increase in the direction of the positive strain (compression). Hence 
the direction of the slow axis will give the direction of the positive strain. 
The retardance will be a projection of the strain in real space to the xy- 
image plane. A quantitative calculation of the strain using the photoelastic 
properties and known crystal orientation [231 El], is °ut of the scope of the 
current proof of concept report. 

3.2.1 Strain imaging of stressed CaF2 

The capability of the MME system to image strain patterns was investigated 
by inducing a direct force onto a crystalline CaF2 prism, which is isotropic 
and transparent from the UV to the infrared. The prism was placed on a 
platform in the image plane of the system, and stress was applied through a 
metal plug, as shown in Figure [4} The upper edge of the prism was opaque, 
such that the interface region between the metal plug and the crystal could 
not be examined. 
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Figure 4: Schematic drawing of the calcium fluoride crystal sample with a 
metal plug causing the observed strain. 



Figure [5] shows Mueller matrix images of the prism without and with 
applied stress. It is evident from the non stressed image that the crystal is 
isotropic, and has no retardance or diattenuation. After the stress is applied, 
the crystal becomes anisotropic, seen by the appearance of the non diagonal 
lower right 3x3 matrix. From the measured Mueller matrix, the magnitude 
of the retardance and the orientation of the slow axis is then calculated. 

The corresponding retardance plots for the two Mueller matrix images in 
Figure [5] are shown in Figure [6} It is evident that when no stress is applied, 
there is no slow axis, and the calculated direction is therefore random, as 
a result of noise. When stress is applied, the direction of the slow axis 
corresponds to the direction of the positive strain. 

This sample could evidently have been studied with any visual MME 
imaging system. However, in case of more rough surfaces or a low band-gap 
crystal, the use of near infra-red light becomes imminent. 

3.2.2 Strain imaging in a multi-crystalline silicon wafer 

The band-gap of Si makes the material transparent for photons with lower 
energy than l.leV, which means that a visible MME can not perform bulk 
strain measurements in transmission. 

In this experiment the sample was initially a 150/xm thick multi-crystalline 
wafer, with a surface of about 4cm x 2.5cm and a trapezoidal shape (see 
Figure [7]). Crystalline silicon is isotropic with a cubic lattice structure. In 
multi-crystalline wafers the grains are through-going and the grain bound- 
aries go straight through the wafer, elliminating the problem of shading 
grains. In order to avoid large losses of intensity and to obtain a good signal 
to noise ratio, the wafer was polished on both sides, down to a thickness of 
approximately 100/im. The wafer was placed in an assembly consisting of 
two bars with slits, keeping the wafer vertically stable and in position (see 
Figure [7| . Stress was then applied to the wafer by buckling [25] , which is a 
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Figure 5: Mueller matrix images of non stressed (a) and stressed (b) calcium 
fluoride crystal. The elements are normalized to my, except from the my 
image showing the intensity image. In (a) the crystal is found isotropic, i.e. 
No strain induced birefringence is observed in the crystal, while in (b) the 
crystal is found with a graded anisotropy due to the non uniform applied 
stress. 
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Figure 6: Calculated retardance from the Mueller matrix image in Figure [5j 
(a) shows the retardance map when no stress is applied, and (b) when stress 
is applied. The lines show the direction of the slow axis, while the color-map 
indicates the retardance in degrees with wavelength of 980nm. 
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Figure 7: The assembly made for the buckling of the wafer. The wafer is 
placed between a fixed and a non fixed bracket, which then can be translated 
by a fine threaded screw. 
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well known method for material strength tests. The strain is induced when 
the bars are forced towards each other such that the wafer bends by the 
control of a fine threaded screw. 

In order to map the whole wafer, the assembly is mounted on a transla- 
tion stage which automatically moves to the next position after each Mueller 
matrix image acquisition. For simplicity and proof of concept, only one of 
the images is considered here, the location of this area on the wafer is indi- 
cated by the red square in Figure [7| 

Figure [8] shows the calculated retardance and the orientation of the slow 
axis of the selected area of the strained wafer. The crystal grains show 
contrast in both the orientation of slow axis and the retardance maps, which 
indicates that both different magnitude and orientation of strain can be 
found. The two areas (1 and 2) indicated by the red squares in Figure [8^a) 
are enlarged and shown in Figure [9] and [TUJ 




Figure 8: Figure (a) shows the calculated retardance map of the area indi- 
cated by the red square in Figure [7j with applied strain, while Figure (b) 
shows the orientation of the slow axis for the area in interest. Area 1 and 
2 indicate the two areas which will be investigated in more detail (Figure [9] 
and [TO]). 



The retardance image of area 1 (Figure |9j) shows a grain boundary where 
the retardance is higher for the grain located in the lower part of the image, 
compared to the upper grain. Various domains of lower and higher retar- 
dance can also be observed inside the two grains. Figure |9^a) shows that at 
the grain boundary there appears to be domains with lower retardance than 
the surrounding grains. It is observed from Figure [9^b), that the orientation 
of the strain is very different in this particular domain. These observations 
may possibly be due to a relaxation of strain at the grain boundary. 

In area 2 (Figure 10) the grain is observed to have a higher retardance 
than the adjecent grains. One may expect that the whole grain experienced 
the same retardance, but in this case a higher retardance is found in the 
lower part of the grain. This observation can possibly be explained by 
the low retardance on the grain boundary around the upper part of the 
grain, where a relaxation may have occurred. The corresponding orientation 
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Figure 9: Figure (a) shows the calculated retardance map of area 1, with 
applied strain. There is good contrast at the grain boundary, also interest- 
ing is the low retardance domain along the grain boundary, which may be 
induced due to a relaxation of strain at the boundary. Figure (b) shows the 
calculated orientation of the slow axis. 



of the slow axis (Figure [lO^b)) shows a different strain orientation than 
the surrounding grains. The details of the structures within the domains 
can without great difficulty be studied in a lab set-up. A production line 
equipment can be evisaged in order to inspect strain over large areas, or 
strain above a particular threshold level. Polarimetric imaging may also 
be useful in several other occasions e.g. to characterize the polarization 
dependent optical properties of optical components, such as lenses, wherever 
these have to be taken into account. 

4 Conclusions 

In conclusion, we have constructed a near infra-red (NIR) Mueller matrix 
imaging ellipsometer based on ferroelectric liquid crystal retarders. The 
system has been demonstrated to be a valuable tool in the characterization 
of strain in NIR transparent solids, and in particular multi-crystalline silicon. 
Through the use of Mueller matrix analysis tools, a detailed analysis of strain 
field domains can be made, which in the future can be quantified and used 
as input to numerical models. 
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Figure 10: Figure (a) shows the calculated retardance map of area 2, with 
applied strain. Figure (b) shows the corresponding orientation of slow axis. 
In this area of the wafer we also see that there is a domain along the grain 
boundary where the retardance is lower than in the two adjacent grains. 
Also the slow axis has a different orientation in this domain. 



imaging system. 
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